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Abstract. Observational data from supernovae type la, baryon acoustic oscillations, gas mass 
fraction in galaxy clusters, and the growth factor are used to reconstruct the the interaction rate 
of the holographic dark energy model recently proposed by Zimdahl and Pavon 1 1 ] in the redshift 
interval < z < 1.8. It shows a reasonable behavior as it increases with expansion from a small or 
, vanishing value in the far past and begins decreasing at recent times. This suggests that the equation 

of state parameter of dark energy does not cross the phantom divide line. 
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INTRODUCTION 

There is a wide-shared conviction among cosmologists that the Universe is nowadays 
^ ■ experiencing a stage of accelerated expansion not compatible with the up to now fa- 

Tij- ! vored Einstein-de Sitter model. This consensus, however, does not extend to the agent 

\ (very frequently called dark energy) behind this acceleration and, at the moment, there 

are many competing candidates [2], the cosmological constant being the leading one. 
However, the latter suffers from two main short-comings: the huge value predicted by 
CX5 ! quantum field theory estimations, and the coincidence problem, i.e., the fact that the 

O ■ energy densities of non-relativistic matter and dark energy are currently of the same or- 

der. This is why many researches are looking for alternative candidates of dark energy. 
\ Among the most recent generic proposals there is a very suggestive one based on the 

^ I holographic principle. Loosely speaking, the latter asserts that the entropy of a system 

is given by the number of degrees of freedom lying on the surface that bounds it, rather 
than in its volume The roots of this principle are to be found in the thermodynamics 
of black holes [4]. Nevertheless, as observed by Cohen et al. ^5*], a system may fulfill the 
holographic principle and, however, include states for which its Schwarzschild radius is 
larger than its size, L. This can be avoided by imposing that the energy of the system 
shouldnot exceed that of black hole of the same size or, equivalently, p < 3c^/ (8;r GL^), 
where is a (non-necessarily constant) parameter. In the cosmological context L, the 
infrared cutoff, is usually taken either as the event horizon radius or the Hubble radius. 
For a quick summary of holographic dark energy see section 3 of Ref. 

The model of Zimdahl et al. [HI] rests in two main assumptions: (z) dark energy 
complies with the holographic principle with L identified as the radius of the Hubble 
horizon, H^, hence = 3c^H^ /{SkG), and (ii) dark energy and dark matter do not 
evolve separately but they interact with one another. As a consequence, the energy 



conservation equations read 



Pm + 3Hp„ = Q, p^ + 3H{l+ w)p, = -Q , (1) 



where w = p^/ px is the equation of state parameter of dark energy which is not con- 
strained to be a constant. Subscripts m and x are for non-relativistic dark matter and dark 
energy, respectively. 

It is to be noted that for spatially flat universes in the absence of interaction, 2 = 0, 
there would be no acceleration [1, 6]. Moreover, Q must be a positive-definite quantity if 
the coincidence problem is to be solved [0] or at least alleviated [[8(], and the second law 
of thermodynamics to be fulfilled |@]. Further if Q were negative, px would have been 
negative in the far past. Besides, it has been forcefully argued that the Layzer-Irvine 
equation [10] when applied to galaxy clusters reveals the existence of the interaction 



UM- To the best of our knowledge, the interaction hypothesis was first introduced, much 



earlier of the discovery of late acceleration, by Wetterich [1121] to reduce the theoretical 
huge value of the cosmological constant, and was first used in the holography context by 
Horvat [13]. As we write, the body of literature on the subject is steadily growing -see 
[Hi [13] and references therein. Most cosmological models implicitly assume that matter 
and dark energy couple gravitationally only. However, unless there exists an underlying 
symmetry that would set Q to zero (such a symmetry is still to be discovered) there is 
no a priori reason to discard the interaction. Ultimately, observation will tell us whether 
the interaction exists. 

Following we will write the interaction as 2 = Tpx, where F is an unknown, 
semi-positive definite, function that gauges the rate at which energy is transferred from 
dark energy to dark matter. Clearly, as long as the nature of both dark ingredients of 
the cosmic substratum remain unknown, F cannot be derived from first principles. The 
alternative is to resort to observational data (in our case, supemovae type la (SN la), 
baryon acoustic oscillations (BAO), gas mass fraction in galaxy clusters and the growth 
factor) to roughly reconstruct it -for details see [15.] . 



RECONSTRUCTING THE RATE 

From Eqs. ([B and the above expressions for px and Q the evolution equation 



r={l+r) 



l+r 



(2) 



for the ratio r = p,,,/ px between the energy densities, is readily obtained. With the help 
of Friedmann equation Q.^ + Q.x + = 1j in terms of the usual density parameters 
= SnGpi/{3H-^) (i = m,x), and Q-k = —k/{a^H'^), where k stands for the spatial 
curvature index of the Friedmann-Robertson- Walker metric, we can write 

r=-2H^q, (3) 



where q = —d/{aH^) denotes the cosmic deceleration parameter. Here, for the holo- 
graphic dark energy we have adopted the expression Px<^ H^. The latter follows from 
choosing the infrared cutoff, L, as the Hubble radius, // ^ 

Likewise, starting from the first of Eqs. ([U) and using Friedmann equation, we get for 
the equation of state parameter 



w{z) = (1+r) 



2H' 
3^" 



3^^ 



(4) 



where z denotes the redshift and a prime indicates derivative with respect to this quantity. 
We fit the Chevallier-Polarsky-Linder parametrization namely, 

w(z) = wo + wi—— , (5) 
1 +z 

where wq is the present value of w{z), and wi a further constant, to current data from 
different observational probes and subsequently use the fitting values for wq and wi to 
reconstruct the dimensionless ratio Y /3H. 

As for the data, we resort to the various SN la observations in recent times. In 
particular we use 60 Essence supemovae [.17.1 . 57 SNLS (Supernova Legacy Survey) 
and 45 nearby supemovae. We have also included the new data release of 30 SNe la 
spotted by the Hubble Space Telescope and classified as the Gold sample by Riess et al. 



[|17ll . The combined data set can be found in Ref. [|18|1 . The total number of data points 
involved is 192. 

Next we add the measurement of the cosmic microwave background (CMB) acoustic 
scale at zbao = 0-35 as observed by the Sloan digital sky survey (SDSS) for the large 
scale structure. This is the BAO peak |[l9|]. 

We also consider the gas mass fraction of galaxy cluster, fgas = Mgas/Mtot, inferred 



from the x-ray observations 112011 . This depends on the angular diameter distance to 

3 /2 

the cluster as fgas = . The number of data points involved is 26. 

Likewise, the two-degree field galaxy redshift survey (2dFGRS) has measured the 
two point correlation function at an effective redshift of Zs = 0.15. This correlation 
function is affected by systematic differences between redshift space and real space 
measurements due to the peculiar velocities of galaxies. Such distortions are expressed 
through the redshift distortion parameter, /3 . Correlation function can be used to measure 
it as /3 = 0.49 ±0.09 at the effective redshift of z = 0.15 of the 2dF survey. This result 
can be combined with linear bias parameter b = 1 .04 ±0.11 obtained from the skewness 
induced in the bispectrum of the 2dFGRS by linear biasing to find the growth factor g at 
z = 0.15, namely ^ = 0.51 ±0.11 



The spatially flat case 



By setting 0.^ = equations Q and (SI) reduce to 



r 

3H 



-ro 



2H' 
3~H 



(6) 



and 



= (l+ro) 



2H' 
3ll~ 



(7) 



respectivel>Q. As usual, indicates the present value of the ratio r. Using these two 
expressions we determine wq and wi from the data and, with them, we reconstruct 
r/3H -see figures [T] and [21 The best fit values, with Icr error bars for the parameters 
when all the data (SN la + BAO + x-rays -i- growth factor) are included, come to 
be: wo = -1.13 ±0.24, wi = 0.66 ± 1.35 (for a,no = 0.25 & Qk = 0, Fig. [B; and 
wo = -0.80 ±0.28, wi = -1.75 ± 1.79 (for a^o = 0.3 & = 0, Fig.©. 

Contrary to what one may think, the fact that r was never large does not bring the 
model of Ref. [1] into conflict with the standard scenario of cosmic structure formation. 
One may believe that at early times the amount of dark matter would have been too 
short to produce gravitational potential wells deep enough for galaxies to condensate. 
However, this is not so; a matter dominated phase is naturally recovered since at high and 
moderate redshifts the interaction is even smaller than at present whence the equation 
of state of the dark energy becomes close to that of non-relativistic matter -see yj for 
details. 

Likewise, it should be noted that dark energy clusters similarly to dark matter when 
the equation of state of the former stays close to that of the latter. To see this more 
clearly we recall the perturbation dynamics of this model studied in [1] by using the per- 
turbed metric d*^ = - ( 1 ± 2 1//) d?^ ± ( 1 - 2 1//) dapdx"dxl^ with yr the scalar metric 
perturbation, and the Bardeen gauge-invariant variable lf22ll 



1 P 



3p+p 



-\j/-H' 



(8) 



The latter represents curvature perturbations on hypersurfaces of constant energy den- 
sity. For the simplest case T — constant, one follows 



Tr 



1 



3r3Hir-T 



(9) 



where the subscript / signals some initial time and a hat indicates perturbation of the 



^ The corresponding equations (namely (6) & (7)) in f\3\ bear an extra factor, (1 This was a typo 
with no further consequences to the calculations. This corrects the typo. 



0.5 



0.5 





FIGURE 1. The dimensionless ratio r/{3H) vs redshift. In the four panels we have fixed n„,o = 0.25 
and £2^. — 0. The soHd Hne is for the mean value and the shaded area indicates the 1 a region. The region 
above the horizontal dashed line can be visited only when the dark energy becomes of phantom type, i.e., 
w < -1. 



corresponding quantity -see Ol] and |15] for details. Accordingly, as far as the the 
equation of state parameter of dark energy w remains close to that of dark matter, w 0, 
both components cluster similarly. 



Non-spatially flat cases 

It is immediately seen that for Q.;^ the ratio r between energy densities is not a 
constant. This introduces a further unknown function in our fitting procedure. Notwith- 
standing, one should not expect a large variation in r in the redshift range (0, 1.8). In 
our subsequent computation, we Taylor expand r around its present value up to the first 
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FIGURE 2. Same as Fig. 1 except that here we have fixed Q.„,q = 0.30 and iljt = 0. 



order term. Accordingly we parameterize it as 



r = ro + ri 



l+z 



(1 



Here ri is a constant which can be related to the present ratio between and by 



This can be used to fix the unknown constant ri for a given Q.i^q and Qjco. Also 



1 f 

r wr 



(1 



m l+r 



(1 



where w is given by Eq. dH). 



With the help of these expressions, the ratio T/SH can be reconstructed from the 
data. The corresponding figure (with Q.„^q = 0.30 and small ^^q) is very similar to Fig. 
2 thereby we do not reproduce it here (see Fig. 3 in Ref. [fisl]'). The best fit values, with 
la error bars for the parameters when all the data (SN la + BAO + x-rays + growth 
factor) are included, come to be: wq = — 0.806±0.29, wi = — 1.74±3.33. We conclude 
that a small spatial curvature, in agreement with the WMAP Syr experiment [.23.1 , has a 
negligible impact on the evolution of the interaction rate. 



DISCUSSION 

Using the observational data (SNIa, BAO, gas mass fraction, and growth factor) available 
in the redshift range < z < 1.8 we reconstructed the interaction term Q of Ref. 01]. 
The interaction rate F (and hence Q) stays positive in the said range. Its general trend 
is to decrease as z increases but it shows no indication of becoming negative at larger 



redshifts. This corroborates that, as previously suggested |!9|,lll|], the energy transfer 
proceeds from dark energy to dark matter, not the other way around. Although phantom 
behavior cannot be excluded at recent and present times it only occurs manifestly either 
for large Q.xo -see Fig.[T]- or when just the supemovae data are considered (top-left panel 
of Figs. [Hand O. When Q.t(0 is somewhat lower (say, 0.7) and BAO and other data are 
included, the mean value of dimensionless interaction rate, r/3H, no longer crosses the 
phantom divide (i.e., the horizontal dashed line). It reaches a maximum near z = to 
subsequently decrease with expansion. This result is rather comforting since holography 
does not seem compatible with phantom energy f24\ . On the other hand, it should be 
noticed that 0.^0 values as high as 0.75 do not appear favored from a combination of 



results from WMAP lyr and weak lensing which yields 0.^0 = 0.70 ±0.3 [|25|] . 

Likewise, a small curvature term -of either sign- is of little consequence. 

At any rate, it is fair to say that the concordance ACDM model (wq = — 1, wi = 0) 
shows compatibility within Icr confidence level with the set of data considered in this 
paper. 
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